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We study the critical temperature T c of SFF trilayers (S is a singlet superconductor, F is a ferromagnetic 
metal), where the long-range triplet superconducting component is generated at noncollinear magnetizations 
of the F layers. We demonstrate that T c can be a nonmonotonic function of the angle a between the mag- 
netizations of the two F layers. The minimum is achieved at an intermediate a, lying between the parallel 
(P, a = 0) and antiparallel (AP, a — n) cases. This implies a possibility of a "triplet" spin-valve effect: at 
temperatures above the minimum Tj 1 but below and T* P , the system is superconducting only in the 
vicinity of the collinear orientations. At certain parameters, we predict a reentrant T c (a) behavior. At the 
same time, considering only the P and AP orientations, we find that both the "standard" (Tjf < T^ p ) and 
"inverse" [Tf > T^ p ) switching effects are possible depending on parameters of the system. 
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In superconducting spin valves with the layer se- 
quence F1/S/F2 the superconducting transition tem- 
perature T c of the system can be controlled by mutual 
alignment of magnetizations Mi^ of the two ferromag- 
netic layers Fl and F2. Therefore, at a temperature T 
fixed inside the range of T c variation, there is an oppor- 
tunity for switching the superconductivity on and off 
by reversing the magnetization direction of the Fl or 
F2 layer. Model calculations have shown that the tran- 
sition temperature p for the antiparallel (Mi M2) 
orientation of the Fl and F2 magnetizations should be 
higher than the transition temperature T p for the op- 
posite case (Mi ft M 2) BSE]- The situation with 
this order of T c 's (i.e., T p < T AP ) is commonly referred 
to as the "standard" switching (see, e.g., [1]), and the 
switching in this case actually occurs at temperatures 
T such that T p < T < T C AP . The basic physical reason 
for the difference AT C = — T c p > is partial com- 
pensation of the pair-breaking ferromagnetic exchange 
field, if the magnetizations of the Fl and F2 layers are 
aligned antiparallel. 

Several experimental groups have published results 
on superconducting spin valves of the F1/S/F2 type 
gl El [El [U [El [HI HI1 H21 H3] - The experimental results 
turned out to be controversial. Some studies of Fl/S /F2 
structures have shown the standard spin-valve effect 
H M M \H IE] with the maximum shift AT C » 41 mK 
reported for the Ni/Nb/Ni trilayer in [? . However, 
some experiments revealed the "inverse" spin-valve ef- 



fect 0Q3J1Q31H3] with T p > T C AP (i.e., AT C < 0). The 
most advanced calculations within the proximity effect 
theory, which take into account the triplet components 
of the superconducting pairing |14j . demonstrate only 
the standard switching [3J 115) with T c monotonically in- 
creasing from the P to AP configuration [3] . Additional 
physical mechanisms like spin imbalance effect [TT] or 
magnetic domain structure [101 112] should be recruited 
to explain the inverse spin-valve effect in the studied 
Fl/S/F2-type structures. 

A bit earlier an unconventional spin-valve-like 
S/F1/F2 structure was theoretically proposed in [IB] 
to control the superconducting T c in the S layer by 
mutual alignment of the magnetizations of the two 
adjacent ferromagnetic layers Fl and F2. The authors 
of [16] argued that T p < T^ p in their system because 
of partial cancelation of the pair-breaking exchange 
fields just within the magnetic F1/F2 subsystem of the 
structure, thus predicting the standard switching as in 
the interleaved F1/S/F2 structure. 

The S/F1/F2 structures are much less investigated 
experimentally [8j 117) , and the experiments indicate the 
standard switching effect [16] with the maximal size of 
about 200 mK. In this Letter we study the critical tem- 
perature of a S/F1/F2 trilayer at arbitrary angle be- 
tween the in-plane magnetizations of the ferromagnetic 
layers (see Fig. 1). We demonstrate that this struc- 
ture allows not only the standard but also inverse spin- 
switching effect. Moreover, we show for the first time 
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Fig. Q] S/F1/F2 trilayer. The S/Fl interface corre- 
sponds to x = 0. The thick arrows in the F layers 
denote the exchange fields h lying in the (y, z) plane. 
The angle between the in-plane exchange fields is a. 

that the minimal critical temperature Tj 1 of the struc- 
ture is achieved at a noncollinear alignment of the mag- 
netizations, when the long-range triplet component of 
the superconducting pairing is generated. Since Tj r is 
lower than both T c p and T^ p , this offers a possibility of 
a "triplet spin-valve effect" never reported before. 

1. Model. We consider the S/F1/F2 structure in 
the dirty limit, which is described by the Usadel equa- 
tions. Near T c , the Usadel equations are linearized and 
contain only the anomalous Green function / [T3J [T5] : 



2 dx 2 



i sgnw 



{f (h<T),/} + Afi<7 =0. (1) 



Here, / is a 4x4 matrix, fj and &i are the Pauli matri- 
ces in the Nambu-Gor'kov and spin spaces, respectively, 
D is the diffusion constant, and ui — irT c (2n + 1) with 
integer n is the Matsubara frequency. The exchange 
field in the middle Fl layer is along the z direction, 
h = (0, 0, h), while the exchange field in the outer F2 
layer is in the yz plane: h = (0, h sin a, h cos a). The 
angle a changes between (parallel configuration, P) 
and 7r (antiparallel configuration, AP). The order pa- 
rameter A is real-valued in the superconducting layer, 
while in the ferromagnetic layers it is zero. In general, 
the diffusion constant D acquires a proper subscript, S 
or F, when Eq. ((TJ is applied to the superconducting 
or ferromagnetic layers, respectively. However, for sim- 
plicity we take them equal in this paper, because this 
assumption does not influence qualitative behavior of 
T c (a). 

The Green function / can be expanded into the fol- 
lowing components: 



/ = h (/boo + h&z + hai) ) 



(2) 



where /q is the singlet component, is the triplet 
with zero projection on the z axis, and f% is the triplet 



with ±1 projections on z (the latter is present only if 
a =/= 0,7r). The singlet component is even in frequency 
(and real- valued) , while the triplet ones are odd (and 
imaginary): /o(-w) = fo(w), h{-w) = -fo(u), and 
f 2 (—uj) — — fi{w), which makes it sufficient to consider 
only positive Matsubara frequencies, u) > 0. 

As we show below, the problem of calculating T c 
can be reduced to an effective set of equations for the 
singlet component in the S layer: the set includes the 
self-consistency equation and the Usadel equation, 



Dd 2 f 



2 dx 2 

with the boundary conditions 
dfo 



UJ>0 v 

w/ + A = 0, 



dx 



= 



dx 



(3) 
(4) 

(5) 



x=0 



Here T c g and £ = \J D/2itT c s are the superconducting 
transition temperature and coherence length for an iso- 
lated S layer, and we assume that the S layer occupies 
the region — ds < x < (see Fig. 1). This is exactly the 
problem for which the multi-mode solution procedure 
(as well as the fundamental-solution method) was de- 
veloped in [TO] and then applied to F1/S/F2 spin valves 
in [3]. We only need to determine the explicit expres- 
sion for W in Eq. ([5]), solving the boundary problem for 
the S/F1/F2 structure. 

2. Solution of the model. To simplify deriva- 
tions, while keeping the essential physics, we consider 
the middle ferromagnetic layer Fl of arbitrary thick- 
ness (0 < x < dp) but the outer ferromagnetic layer F2 
being semi-infinite {dp < x < oo). The Usadel equation 
(HJ) generates the following characteristic wave vectors: 




kh = 



Uk\. 



(6) 



Only k u appears in the solution for the S layer, while 
the F-layers' solutions are described by k u , kh, and k^- 
Since the exchange energy is usually larger than the su- 
perconducting energy scale, h T c , the k u mode in the 
ferromagnetic layers (arising at noncollinear magnetiza- 
tions) represents the long-range triplet component [T4"] . 
which plays the key role in the present study. 
In the S layer the solution of Eq. ([1]) is: 





cosh(fc w (x + ds)) 
cosh (k u ds) 



(7) 
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The singlet component fo(%) in the S layer cannot be 
written explicitly, since it is self-consistently related to 
the (unknown) order parameter A(x) by Eqs. (|3"] ) -([S |) . 
Our strategy now is to obtain the effective boundary 
conditions ([5]) for fo(x), eliminating all other compo- 
nents in the three layers. 

In the middle Fl layer the solution of Eq. (jTJ) reads: 



ffo(x) 

\f2(x) 

+C 2 



+s 2 



/0> 

I cosh (k u x) + Si 

U 



I sinh(fc w a:) + 



cosh 



smb. ( kh 




cosh ( kfrX 



sinh f kl 



(8) 



Finally, the solution in the semi-infinite, outer F2 
layer is built only from descending modes: 



exp (— k u (x - (If)) + 




+ E 3 



cos a | exp 
, sin a 



-kh(x - d F )) + 



-kfrix — dp) 



(9) 



We will use the simplest, perfect-transparency 
boundary conditions at the S/Fl and F1/F2 interfaces 
(the case 7 = 1 and jb — in the notations of 20 ): 



f, 



* I left 



fi\ 



right ' 



dx 



loft 



dx 



right 



Altogether there are 12 boundary conditions at the two 
interfaces (S/Fl and F1/F2). We are mainly interested 
in one of them, determining the derivative of the singlet 
component on the S side of the S/Fl interface (x = 0): 



dfo 
dx 



2Re(fc h 5 2 ) 



(11) 



The remaining 11 boundary conditions form a system 
of 11 linear equations for 11 coefficients entering Eqs. 
©-(E]). The solution of this system is nonzero due to 
/o(0) coming from Eq. ([7]) and entering the "right-hand 




2 kfo dp 

Fig.d Dependence of W(0) -W(ir), Eq. |[T2]). on the 
thickness c?f of the Fl layer. Positive values of this 
oscillating function correspond to stronger suppression 
of superconductivity at the P alignment (the standard 
switching effect), while negative values correspond to 
stronger suppression of superconductivity at the AP 
alignment (the inverse switching effect). 

side" of the system. Finding the S2 coefficient [which is 
proportional to /o(0)], we substitute it into Eq. pT|) and 
thus explicitly find W entering the effective boundary 
conditions (JSJ. 

3. Analysis of the solution. After reducing the 
problem to Eqs. ([3])-([5]), all the information about the 
two F layers is contained in the single real-valued func- 
tion W. This function makes fo(x) bend at the S/Fl 
interface, hence the larger W, the stronger T c is sup- 
pressed. 

The explicit expression for W(a) is very cumber- 
some and we do not write it here. However, certain 
analytical development (as well as complete numerical 
analysis) is possible. For the analytical consideration, 
we make an additional assumption of T c -C h, which 
implies k u <C kh- For the collinear cases (a = and 
a = 7r) we then find W(0) = 2fc/j£ and 



(10) W(0) - W(ir) - 2k h ( 



-2k h d f 



y/2sm(2k h d F + 7r/4) 

sinh^fchdp) + cos(2/cftd_p) 



(12) 

which oscillates as a function of dp, changing its sign 
(see Fig. 2). Thus, as a result of interference in the 
middle Fl layer, we can either have the standard spin- 
switching effect with T P < X^ p [when the pair-breaking 
at the P alignment is stronger than at the AP align- 
ment of magnetizations, i.e., at W(0) — W(tt) > as in 
the range Ik^dp < in/A in Fig. 2] or the inverse spin- 
switching effect with T c p > T C AP [at W{0)-W(ir) < as 
in the range 37r/4 < 2khdp < 77r/4 in Fig. 2]. Note that 
the amplitude of the inverse effect is notably smaller 
compared with the standard one. The analytical calcu- 
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Fig. [3] Critical temperature T c vs. the misalignment 
angle a for various thicknesses of the Fl layer. We took 
h/nT c s = 6.8; all other parameters are shown in the 
figure. In the cases d F — and d F = oo, which are 
physically equivalent (curve "1"), T c does not depend 
on a. Curves "2" and "4" correspond to the standard 
and inverse switching effects, respectively. Curve "3" 
demonstrates the triplet spin- valve effect. The coher- 
ence lengths £s and £f were taken equal (denoted by £ 
in the text) in order to present our main results in the 
simplest possible case. At a = all the curves coincide, 
since in this case the F part of the system is physically 
equivalent to a single half-infinite F layer. 



lation of the second derivatives of W(a) at a — and it 
(the first ones are zero) shows that under the above as- 
sumption, both the collinear alignments represent local 
minima of W{a). This means that T c (a) decreases as 
the configuration deviates from the P or AP alignment. 
Therefore, T c (a) is nonmonotonic, and the minimal T c 
must be achieved at some noncollinear configuration of 
magnetizations at a ^ 0, n. 

The analytical results obtained at k u <C kh are il- 
lustrated and extended by numerical calculations at ar- 
bitrary relation between k u and kh- Figure 3 shows 
dependence of the transition temperature T c on the an- 
gle a between the magnetizations. We see that at small 
thicknesses dp of the middle ferromagnetic layer Fl, the 
switching effect is standard, while at larger dp the ef- 
fect is inverse (T c p > T^ p ). Moreover, when the Fl layer 
thickness is around a half of the coherence length £, the 
minimal critical temperature Tj 1 at noncollinear orien- 
tations is significantly lower that both T p and T^ p - 
this case corresponds to the triplet spin- valve effect. 
Note that depending on the parameters of the system, 
the minimum of T c (a), predicted analytically, can shift 
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Fig-HJ T c (a) at various d F and ds. Curve "1" coincides 
with curve "3" in Fig. 3. Curves "2" and "4" demon- 
strate the reentrant behavior, in which case the triplet 
spin-valve effect takes place even at T = 0. 



to a close vicinity of either a = or a = tt, becoming 
shallow and indistinguishable. 

Figure 4 demonstrates the possibility of reentrant 
T c (a) dependence. In this situation the triplet spin- 
valve effect takes place even at T = 0. 

4. Discussion. The physical interpretation of the 
triplet spin-valve effect can be given as follows: at the 
collinear configurations, both the singlet component fo 
and the zero-projection triplet component fa of the 
pairing function are short-ranged (with the character- 
istic penetration depth of the order of k^ 1 ), so that at 

1 ^ the middle Fl layer plays a role of a shield 
separating the S layer from the ferromagnetic half-space 
F2. When the angle between magnetizations declines 
from the collinear configurations, the long-range triplet 
component f% of the pairing function is generated [T3] . 
Then, the S layer becomes effectively coupled by this 
long-range triplet component to the semi-infinite fer- 
romagnetic F2 layer. The pair-breaking in the S layer 
enhances, giving rise to more effective suppression of su- 
perconducting T c . In other words, we can say that the 
T c suppression is due to "leakage" of Cooper pairs into 
the ferromagnetic part. In this language, the generation 
of the long-range triplet component opens up an addi- 
tional channel for this "leakage" , hence T c is suppressed 
stronger. 

In order to supply the qualitative picture by more 
quantitative details, we can find the amplitudes of dif- 
ferent components at the S/Fl interface in the limit of 
ku) <C kh and large dp. This can be done analytically 
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from the boundary conditions which produce the linear 
system of equations for the coefficients entering Eqs. ([7])- 
(|9]). We find that in the limit of cLf k~ x , k^ 1 , the am- 
plitudes of the long-range triplet components near the 
S/Fl interface (which are given by C±, Si, and B) are 
suppressed by the factor e ~~ ku ' dF ~ khdF which has a clear 
physical interpretation. The long-range components are 
generated from the short-range ones at the F1/F2 inter- 
face (i.e., at x = dp), where electrons "feel" inhomoge- 
neous magnetization. Therefore, the long-range contri- 
bution at the S/Fl interface is obtained as a result of 
a "wave" that goes from the S/Fl interface as a short- 
range component with the wave vector kh and returns 
after reflection at the F1/F2 interface as a long-range 
component with the wave vector k u . At the same time, 
the self-consistency equation (j3|) that determines T c , 
contains only the singlet short-range component. There- 
fore, the influence of the long-range components on T c 
is indirect: the long-range components influence T c only 
through their influence on the singlet component. We 
find that while the difference between W (that encodes 
the information about the suppression of T c ) for the 
P and AP cases is suppressed as e~ 2khdp (the short- 
range components go from the S/Fl to F1/F2 interface 
and back), the changes in W due to noncollinear mag- 
netizations contain the same exponential, e ~ 2khdF . Of 
course, the influence of the long-range triplet compo- 
nents is contained in prcfactors but no long-ranged ex- 
ponential (with k u instead of kh) appears in W, because 
W still originates from the short-range components. 

In conclusion, we have considered a mesoscopic 
S/F1/F2 structure composed of a superconducting layer 
S, a ferromagnetic layer of arbitrary thickness Fl, and 
a ferromagnetic half-space F2. We have demonstrated 
that the structure exhibits different relations between 
the critical temperatures in the parallel and antiparal- 
lel configuration: both the standard (T c p < T^ p ) and 
inverse (T p > T C AP ) switching can be realized depend- 
ing on the system's parameters. At the same time, our 
main result is that Tj 1 at noncollinear magnetizations is 
lower than both TjF and T^ p , which makes this system 
a triplet spin valve. Possible experimental observation 
of a nonmonotonic (like curve "3" in Fig. 3 or curves "1" 
and "3" in Fig. 4) or even reentrant (like curves "2" and 
"4" in Fig. 4) behavior of T c (a) could be a signature of 
existence of the long-range triplet superconducting cor- 
relations [14] in SF hybrid structures. 
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